Odd Triplet Pairing in clean Superconductor/Ferromagnet heterostructures 
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We study triplet pairing correlations in clean Ferromagnet (F)/Superconductor (S) nanojunctions, 
via fully self consistent solution of the Bogoliubov-de Gennes equations. We consider FSF trilayers, 
with S being an s-wave superconductor, and an arbitrary angle a between the magnetizations of 
the two F layers. We find that contrary to some previous expectations, triplet correlations, odd in 
time, are induced in both the S and F layers in the clean limit. We investigate their behavior as a 
function of time, position, and a. The triplet amplitudes are largest at times on the order of the 
inverse "Debye" frequency, and at that time scale they are long ranged in both S and F. The zero 
temperature condensation energy is found to be lowest when the magnetizations are antiparallel. 
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The proximity effects in superconductor/ferromagnet 
(SF) heterostructures lead to the coexistence of ferromag- 
netic and superconducting ordering and to novel trans- 
port phenomena P, Interesting effects that arise from 
the interplay between these orderings have potential tech- 
nological applications in fields such as spintronicsQ- For 
example, the relative orientation of the magnetizations in 
the F layers in FSF trilayers can have a strong influence 
on the conductivity 0, H, H, 0, HI > making them good spin 
valve candidates. Such trilayers were first proposed B for 
insulating F layers and later for metallic! 10. 1 1 XI] ones. 

This interplay also results in fundamental new physics. 
An outstanding example is the existence of "odd" triplet 
superconductivity. This is an s-wave pairing triplet state 
that is even in momentum, and therefore not destroyed 
by nonmagnetic impurities, but with the triplet corre- 
lations being odd in frequency, so that the equal time 
triplet amplitudes vanish as required by the Pauli prin- 
ciple. This exotic pairing state with total spin one was 
proposed long ago [l2| as a possible state in superfluid 
3 He. Although this type of pairing does not occur there, 
it is possible in certain FSF systems (ll. 13. Il3. 14 1 with or- 
dinary singlet pairing in S. This arrangement can induce, 
via proximity effects, triplet correlations with m = and 
to = ±1 projections of the total spin. If the magnetiza- 
tion orientations in both F layers are unidirectional and 
along the quantization axis, symmetry arguments show 
that only the m = projection along that axis can exist. 

Odd triplet pairing in F/S structures has been studied 
in the dirty limit through linearized Usadel-type quasi- 
classical equations 0, O, [3, EH- In this case, it was 
found that to = triplet pairs always exist. They are 
suppressed in F over short length scales, just as the sin- 
glet pairs. The m = ±1 components, for which the ex- 
change field is not pair-breaking, can be long ranged, and 
were found to exist for nonhomogeneous magnetization. 
For FSF trilayersjl, [3, 0], the quasiclassical methods 
predict that the structure contains a superposition of all 




FIG. 1: Schematic of FSF junction. The left ferromagnetic 
layer Fi has a magnetization oriented at an angle —a/2 in 
the x — z plane, while the other ferromagnet, F2, has a mag- 
netization orientation at an angle a/2 in the x — z plane. 



three spin triplet projections except when the magneti- 
zations of the F layers are collinear, in which case the 
to = ±1 components along the magnetization axis van- 
ish. It is noted in Ref. [l[ that the existence of such effects 
in the clean limit has not been established and may be 
doubted. This we remedy in the present work, where 
we establish that, contrary to the doubts voiced there, 
induced, long-ranged, odd triplet pairing does occur in 
clean FSF structures. 

Experimental results that may argue for the existence 
of long range triplet pairing of superconductors through 
a ferromagnet have been obtained in superlattices[l8( 
with ferromagnetic spacers, and in two superconduc- 
tors coupling through a single ferromagnet [l9l l2fij ]. 
Measurements [l9| on a SQUID, in which a phase change 
of 7r in the order parameter is found after inversion, in- 
dicate an odd-parity state. Very recently, a Josephson 
current through a strong ferromagnet was observed, in- 
dicating the existence of a spin triplet state [l^ induced 
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by NbTiN, an s-wave superconductor. 

In this paper, we study the induced odd triplet 
superconductivity in FSF trilayers in the clean limit 
through a fully self-consistent solution of the microscopic 
Bogoliubov-de Gennes (BdG) equations. We consider ar- 
bitrary relative orientation of the magnetic moments in 
the two F layers. We find that there are indeed induced 
odd triplet correlations which can include both m = 
and m = ±1 projections. We directly study their time 



dependence and we find that they are largest for times of 
order of the inverse cutoff "Debye" frequency. The corre- 
lations are, at these time scales, long ranged in both the 
S and F regions. We also find that the condensation en- 
ergy depends on the relative orientation of the F layers, 
being a minimum when they are antiparallcl. 

To find the triplet correlations arising from the non- 
trivial spin structure in our FSF system, we use the BdG 
equations with the BCS Hamiltonian, H e g: 
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where A(r) is the pair potential, to be determined self- 
consistently, ipg , ipg are the creation and annihilation op- 
erators with spin S, Ep is the Fermi energy, and a are 
the Pauli matrices. We describe the magnetism of the F 
layers by an effective exchange field h(r) that vanishes in 
the S layer. We will consider the geometry depicted in 
Fig. [1] with the y axis normal to the layers and h(r) in 
the x — z plane (which is infinite in extent) forming an 
angle ±a/2 with the z axis in each F layer. 

Next, we expand the field operators in terms of a Bo- 
goliubov transformation which we write as: 



ips(r) = 22 (u«5(r)7n + Vsv n s(r)% 
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where 77,5 = 1(— 1) for spin down (up), u n s and v n s are the 
quasiparticle and quasihole amplitudes. This transforma- 
tion diagonalizes H e g: [H e s,7n] = ~e n 7n, [HeS,lt] = 
e n ln- By taking the commutator [if)$ (r) , H e s ] , and with 
h(r) in the x — z plane as explained above, we have the 
following: 

hh(r),Weff] = (n c - M^t(r) - K1>i(r) + A(r)^{(r), 

(2a) 

[^(r),74ff] = (n c + h^r) - h x ^(r) - A(r)^{(r). 

(2b) 

Inserting (TTJ into (2]| and introducing a set p of Pauli- like 
matrices in particle-hole space, yields the spin-dependent 
BdG equations: 

[pz ® (ft i ~ h * a z) + (My)Px - h x l) ® a x ] $„ = e n $„, 

(3) 

where $„ = {u n ^(y), u nl (y), w„ T (y), v nl (y)) T and H = 
—dy/(2m) + e± — Ep. Here e± is the transverse kinetic 
energy and a factor of e lk± r has been suppressed. In 
deriving Eq. |(3} care has been taken to consistently use 
the phase conventions in Eq. (JTJ) . To find the quasiparti- 
cle amplitudes along a different quantization axis in the 



x—z plane, one performs a spin rotation: $„ — > U(a')$ n , 
where U(a') = cos(a'/2)l ® 1 — zsin(a'/2)p z ® a y . 

When the magnetizations of the F layers are collinear, 
one can take h x — 0. For the general case shown in 
Fig. [U one has in the Fi layer, h x — ho sin(— a/2) and 
h z = ho cos(— a/2), where ho is the magnitude of h, 
while in F2, h x = ho sin(a/2), and h z — /ioCOs(a/2). 
With an appropriate choice of basis, Eqs. ((3| are cast 
into a matrix eigenvalue system that is solved itera- 
tively with the self consistency condition, A(y) = g(y)fz 
{h - \ [(^WV-iW) - ^(r)V T (r))]). In the F layers 
we have g{y) = 0, while in S, g(y) = g, g being the usual 
BCS singlet coupling constant there. Through Eqs. fl}, 
the self-consistency condition becomes a sum over states 
restricted by the factor g to within wd from the Fermi 
surface. Iteration is performed until self-consistency is 
reached. The numerical process is the same that was used 
in previous work 24 , 25[ , with now the h x term requiring 
larger four-component matrices to be diagonalizcd. 

We now define the following time dependent triplet am- 
plitude functions in terms of the field operators, 

/o(r,t) =\ [<V T (r,*)^(r,0)> + ^;(r,*)^ T (r,0)>] , (4a) 

h(r,t) =i[^ T (r,i)^(r,0)>-(^(r,i)^(r,0)>], (4b) 

which, as required by the Pauli principle for these s-wave 
amplitudes, vanish at t — 0, as we shall verify. Making 
use of Eq. |T]) and the commutators, one can derive and 
formally integrate the Heisenberg equation of the motion 
for the operators and obtain: 

fo(V,t) =^X^ Un T(y) U "-l(y) - U nl(y)Vnl(y)}(n(t), (5a) 

71 

fi(y,t) =- ]^22[un]{y)vn]{y) +u n i(y)v ni (y)]Q n {t), 

n 

(5b) 
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FIG. 2: (Color online) The real part, fo, of the triplet ampli- 
tude fo, for a FSF trilayer at 7 different times. We normalize 
fo by the singlet bulk pair amplitude, Ao/g. The coordinate 
y is scaled by the Fermi wavevector, Y = &f2/, and time by 
the Debye frequency, r = wot. At r = 0, /o = as required 
by the Pauli principle. The interface is marked by the verti- 
cal dashed line, with an F region to the left and the S to the 
right. Half of the S region and part of the left F layer are 
shown. The inset shows the maximum value of fo versus r. 

where Cn(t) = cos(e„i) — i sin(e„t) tanh(e„/2T). 

The amplitudes in Eqs. (O contain all information on 
the space and time dependence of induced triplet correla- 
tions throughout the FSF structure. The summations in 
Eqs. ([5]) are over the entire self-consistent spectrum, en- 
suring that fo and f\ vanish identically at t = and thus 
obey the exclusion principle. Using a non self consistent 
A(y) leads to violations of this condition, particularly 
near the interface where proximity effects are most pro- 
nounced. Geometrically, the indirect coupling between 
magnets is stronger with fairly thin S layers and rela- 
tively thick F layers. We thus have chosen ds — (3/2)£ 
and dpi — dp 2 — £o, with the BCS correlation length 
£o = lOOfc^ 1 . We consider the low T limit and take 
lue> = O.OAEp. The magnetic exchange is parametrized 
via / = ho/Ep- Results shown are for J = 0.5 (unless 
otherwise noted) and the magnetization orientation an- 
gle, a, is swept over the range < a < -k. No triplet 
amplitudes arise in the absence of magnetism (1 = 0). 

For the time scales considered here, the imaginary 
parts of fo(y,t) and fi(y,t) at t ^ are considerably 
smaller than their real parts, and thus we focus on the 
latter, which we denote by fo(y, t) and fi{y, t). In Fig. [21 
the spatial dependence of fo is shown for parallel mag- 
netization directions (a = 0) at several times r = u}p>t. 
The spatial range shown includes part of the F\ layer 
(to the left of the dashed line) and half of the S layer 
(to the right). At finite r, the maximum occurs in the 
ferromagnet close to the interface, after which fo under- 
goes damped oscillations with the usual spatial length 



FIG. 3: (Color online) Spatial and angular dependence of /i, 
at r = 4 ~ t c and several a. Normalizations and ranges are 
as in Fig. [2] Inset: maxima of fo and fi in Fi versus a. 

scale £/ w (kp-^ — kpi) -1 w kp /I. The height of the 
main peak first increases with time, but drops off after 
a characteristic time, r c « 4, as seen in the inset, which 
depicts the maximum value of fo as a function of r. As 
r increases beyond r c , the modulating fo in F develops 
more complicated atomic scale interference patterns and 
becomes considerably longer ranged. In S, we see imme- 
diately that fo is also larger near the interface. Since 
the triplet amplitudes vanish at r = 0, short time scales 
exhibit correspondingly short triplet penetration. The 
figure shows, however, that the value of fo in S is sub- 
stantial for t > t c , extending over length scales on the 
order of £o without appreciable decay. In contrast, the 
usual singlet correlations were found to monotonically 
drop off from their r = value over r scales of order 
unity. 

In the main plot of Fig. [3] we examine the spatial de- 
pendence of the real part of the m = ±1 triplet ampli- 
tude, f\. Normalizations and spatial ranges are as in 
Fig. [2] but now the time is fixed at t = 4 w r c , and 
five equally spaced magnetization orientations are con- 
sidered. At a = 0, /i vanishes identically at all t, as 
expected. For nonzero a, correlations in all triplet chan- 
nels are present. As was found for fo, the plot clearly 
shows that f\ is largest near the interface, in the F re- 
gion. Our geometry and conventions imply (see Fig. Q]) 
that the magnetization has opposite ^-components in the 
Fi and F2 regions. The f\ triplet pair amplitude profile 
is thus antisymmetric about the origin, in contrast to the 
symmetric fo, implying the existence of one node in the 
superconductor. Nevertheless, the penetration of the f\ 
correlations in S can be long ranged. We find that f\ and 
fo oscillate in phase and with the same wavelength, re- 
gardless of a. The inset illustrates the maximum attained 
values of fo and f\ in Fi as a varies. It shows that for 
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FIG. 4: (Color online) The T = condensation energy, Ai?o, 
normalized by /V(0)Aq (A^O) is the usual density of states), 
vs. the angle a for two values of I. When the two magne- 
tizations are antiparallel (a = n) AEo is lowest. The inset 
shows the ordinary (singlet) pair potential averaged over the 
S region, normalized to the bulk Aq. 



a and are larger at I = 1 than at J = 0.5. The half- 
metallic case of J = 1 illustrates that by having a single 
spin band populated at the Fermi surface, Andreev reflec- 
tion is suppressed, in effect keeping the superconductivity 
more contained within S. 

Thus, we have shown that in clean FSF trilayers in- 
duced odd triplet correlations, with m = and m = ±1 
projections of the total spin, exist. We have used a mi- 
croscopic self-consistent method to study the time and 
angular dependence of these triplet correlations. The 
correlations in all 3 triplet channels were found, at times 
t = uj£>t > t c , where r c w 4, to be long ranged in both 
the F and S regions. Finally, study of the condensation 
energy revealed that the ground state energy is always 
lowest for antiparallel magnetizations. 

This project was supported in part by a grant of HPC 
resources from the ARSC at the University of Alaska 
Fairbanks (part of the DoD HPCM program) and by the 
University of Minnesota Graduate School. 
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a broad range of a, a < 37r/4, the maximum of fo varies 
relatively little, after which it drops off rapidly to zero at 
a = 7r. This is to be expected as the anti-parallel orienta- 
tion corresponds to the case in which the magnetization 
is in the x direction, which is perpendicular to the axis 
of quantization (see Fig. [I]). The rise in the maximum 
of fi is monotonic, cresting at a = n, consistent with 
the main plot. At this angle the triplet correlations ex- 
tend considerably into the superconductor. At a = n/2 
the maxima coincide since the two triplet components 
are then identical throughout the whole space because 
the magnetization vectors have equal projections on the 
x and z axes. At a = ir both magnetizations are normal 
to the axis of quantization z (see Fig.[TJ. By making use 
of the rotation matrix U (see below Eq. [3|) one can verify 
that the m = ±1 components with respect to the axis x 
along the magnetizations are zero. 

We next consider the condensation energy, AE , cal- 
culated by subtracting the zero temperature supercon- 
ducting and normal state free energies. The calculation 
uses the self consistent spectra and A(y), and methods 
explained elsewhere 25j, [26j . In the main plot of Fig. [4] 
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we show AEo (normalized at twice its bulk S value) at 
two different values of /. The condensation energy results 
clearly demonstrate that the antiparallel state (a — 7r) is 
in general the lowest energy ground state. These results 
are consistent with previous studies @] of FSF structures 
with parallel and antiparallel magnetizations. The inset 
contains the magnitude of the spatially averaged pair po- 
tential, normalized by Ao, at the same values of /. The 
inset correlates with the main plot, as it shows that the 
singlet superconducting correlations in S increase with 
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